Ultrafast dynamics of surface electromagnetic waves in nanohole array on metallic film 
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We study the ultrafast dynamics of surface electromagnetic waves photogenerated on aluminum 
film perforated with subwavelength holes array by means of transient photomodulation with ~ 100 fs 
time resolution. We observed a pronounced blueshift of the resonant transmission band that reveals 
the important role of plasma attenuation in the dynamics and that is inconsistent with plasmon- 
polariton mechanism of extraordinary transmission. The transient photomodulation spectra were 
successfully modeled within the Boltzmann equation approach for the electron-phonon relaxation 
dynamics, involving non-equilibrium hot electrons and quasi-equilibrium phonons. 



I. INTRODUCTION 



On a smooth metal-dielectric interface light does not 
couple to elementary electronic excitations of the metal 
surfaces because of energy and momentum conserva- 
tion. But in a metal film that is perforated with a two- 
dimensional (2D) periodic array of holes forming a 2D 
metallo-dielectric photonic crystal, the periodicity pro- 
motes zone folding that results in the formation of band 
structure. The grating makes it possible for light to cou- 
ple directly to surface electromagnetic waves by trans- 
fcring them the reciprocal lattice momentum g = 2n/a, 
where a is the grating period^. Indeed, it was recently 
found 2 -! 3 -^ that the optical transmission through hole ar- 
rays fabricated on optically thick metallic films is en- 
hanced by orders of magnitude at resonances; this phe- 
nomenon has been dubbed " extraordinary transmission" 
(EOT). The maxima in the zero-order transmission spec- 
trum is alternated by deep minima caused by Wood's 
anomalies formed due to the diffraction grating of the 
corrugated him. The optical properties of various sub- 
wavelength hole array films have been widely studied 
and different models have been proposed to explain the 
underlying mechanism of the anomalous transmission. 
While earlier models atributed the enhancement to exci- 
tation of film's surface plasmon-polaritons (SPP ffi^a^ , 
more recent works adrresed the role resonant cavity 
modea 7 ^^, dynamical diffraction resonance s 10 ' 11 , and in- 
terference of diffracted evanescent wavesi 2 -. While the 
precise mechanism of extraordinary transmission is still 
under discussion, the transient dynamics of surface elec- 
tromagnetic modes might shed more light on the origins 
of this phenomenon. 

Here we investigate the ultrafast dynamics of surface 
electromagnetic modes measured by using the pump- 
probe correlation spectroscopy of the transmission bands 
in 2D subwavelength hole array on aluminum (Al) thin 
film. We found significant difference with the dynam- 
ics observed previously in smooth metal films (mostly 
gold and silver) in both the time-resolved and spectrally- 
resolved signals. Following an instantaneous rise at 
the onset of the impinging pulse, the transient differen- 
tial transmission exhibits a fast rise on subpicosecond 



timescale followed by a plateau for time delays of 1-3 ps, 
and subsequent slow decay with characteristic time of 
~ 40 ps. The short-time dynamics is accompanied by a 
blueshift of the main transmission band in the transient 
spectra that persists for time delays as long as 100 ps. 
The observed dynamics can be explained by a fast en- 
ergy transfer in the Al film from the electron gas to the 
lattice, with subsequent cooling of the Al film by heat 
transfer to the glass substrate. The fast lattice tempera- 
ture rise is caused by strong electron-phonon interaction 
in Al, which makes the electron-lattice energy transfer 
rate comparable to the rate of nonequilibrium electrons 
thermalization via electron-electron interactions 13 . Im- 
portantly, the blueshift is incosistent with the SPP mech- 
anism and points to the interference as the source of the 
extraordinary transmission. 



II. EXPERIMENT 

The plasmonic lattice sample was a 70 nm thick Al 
film, 5x5 mm 2 in area deposited on a glass substrate. 
The metal film was perforated with circular holes having 
diameter D = 150 nm in a square hole array with lattice 
constant a = 300 nm resulting in a fractional aperture 
area of 25%. The ultrafast laser system used for measur- 
ing the transient PM spectrum was a Ti:sapphire regen- 
erative amplifier with pulses of 100 fs duration at photon 
energies of 1.55 eV, with 400 energy per pulse at a 
repetition rate of 1 kHz. The second harmonic of the 
fundamental pulses at 390 nm (3.1 eV) were used as the 
pump beam. The probe beam was a white light super- 
continuum generated in a 1-mm thick sapphire plate that 
covers the spectral range from 1.6 to 2.8 eV; the white 
light super-continuum was also used for measuring the 
unperturbed transmission spectrum Tq(u>). The pump 
and probe beams were directed onto the perforated Al 
film from the air side to maximize coupling of light to the 
SEW. The transient PM spectrum was obtained from the 
photoinduced change (AT) in Tq using a phase-sensitive 
technique with a resolution AT/Tq ~ 10~ 4 that corre- 
sponds to a photoexcited electron density in the Al metal 
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of <~ 10 cm 3 /pulse. 
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FIG. 1: The measured (solid line) and calculated (dashed 
line) normal incidence optical transmission, To through an Al 
plasmonic lattice comprised of a square array of nanoholes 
with lattice constant 300 nm and hole diameter 150 nm. The 
inset shows the sample AFM image. Numerical fit is based 
on the effective medium model for linear transmission in Sec. 

MM 



Figure [TJ shows the EOT spectrum of the Al plas- 
monic lattice at normal incidence. Several EOT bands 
and transmission anti-resonances (AR) are observed. We 
note that Tq(lu) at resonance frequencies substantially 
surpasses the fractional aperture area of 25%. The low- 
est frequency EOT band peaks at ~ 2.1 eV; whereas 
the lowest frequency AR feature that is associated with 
Wood's anomaly at the Al/glass interface with wavevec- 
tor £;( 01 ) = ge l^ 2 (where e g i — 1.56 and g = 2n/a ) is 
in the form of a dip at 2.6 eV. The AR's at higher fre- 
quencies are higher-order Wood's anomalies that corre- 
spond to wavevectors A;( m ") = g (m 2 + n 2 ) 1 ^ 2 e d (m, 
n are integers) as follows: 3.7 eV (fc( n ) = ge gi 1 ^ 2 2 1 / 2 ) 
that corresponds to the (1,1) branch at the Al/glass in- 
terface; and 4.0 eV (fc( 01 ' = g) that corresponds to the 
(0,1) branch at the Al/air interface. It is noteworthy that 
the lowest frequency EOT resonance (at ~ 2.1 eV) oc- 
curs at considerably lower (by ~ 0.6 eV) frequency than 
that predicted by the lowest-order SPP with wavevector 
satisfying the Bragg condition 2 

fci™")=5[(m 2 + n 2 )( £ - 1 + £d - 1 )] 1/2 . (1) 

The transient PM spectra of the perforated Al film 
are shown in Figs. [2] - HI The time-evolution of pump- 
probe signal at 650 nm, corresponding to probe photon 
energy of 1.9 eV, reveals two component of the dynam- 
ics. Firstly, there is a fast rise of pump-probe signal on 
timescale < 1 ps, followed by a plateau at 1-3 ps (see Fig. 

Secondly, after about 3 ps the saturated PM decays 
with characteristic time of ~ 40 ps (see Fig. [3]) due the 
heat transfer to the glass substrate. Figure @[a) shows 
the transient PM spectral change from 1.7 to 2.4 eV cor- 
responding roughly to the main transmission band spec- 
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FIG. 2: Short-time evolution of differential transmission at 
1.9 eV is shown together with numerical fit (solid line) based 
on solution of the Boltzman equation in Sec. MI Bl 



-5 - . 



1— 




1— 


-10 - 


< 




CO 




o 


-15 - 






-20 - 




-25 - 



Time delay (ps) 

FIG. 3: Long-time evolution of differential transmission at 1.9 
eV due to heat transfer to the glass matrix is shown together 
with numerical fit (solid line) based on two-temperature 
model in Sec. IHTBI 

tral width. A transient blueshift is apparent that persists 
for times up to 100 ps, i e. much longer that electron gas 
thermalization time, pointing to a prominent role of e-ph 
interactions in the dynamics. However this shift can- 
not be well reproduced by a simple spectral derivative 
[see Fig. HJb)] indicating that a perturbation approach 
cannot properly describe the observed changes in optical 
spectra for our excitation intensities. At the same time, 
although the data becomes more noisy at higher energies 
as one approaches to the Wood anomaly, a clear nonlin- 
earity in the frequency dependence emerges. The energy 
dependence of dT/T is determined by transient change 
in the mode dispersion that reflects the periodicity of the 
structure. 

The short-time dynamics exhibts significant differences 
with that observed in gold and silver filrn a 14 ' 15 . In these 
studies, the picosecond dynamics consisted of two well- 
separated in time stages: thermalization of nonequilib- 
rium electrons due to electron-electron (e-e) scatering 
at 200-500 fs accompanied by the signal rise, followed 
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by the decay on a 1-2 ps timescale due to the cooling 
of hot electron gas to the metal lattice. The dominant 
factor that determines the transient dynamics in noble 
metals is the change in the interband succeptibilily of 
the highly polarizable ci-band electrons that is very sen- 
sitive to the smearing of the Fermi edge due to the rise 
in electron temperature T e . In contrast, in simple met- 
als such as Al, the transient optical spectra are gov- 
erned by the changes in the intraband dielectric func- 
tion £ m (u>) = 1 — lo 2 p /uj(u) + ij), where oj p is the bulk 
plasmon energy and 7 is the optical relaxation rate^. 
The absorbance of the film and hence the transient dy- 
namics is related to the deviation of 7 from its intrin- 
sic value 70 that mainly originates from the change in 
the electron-phonon (e-ph) scattering, while the contri- 
butions from e-e, and impurity scatering are relatively 
weakii. Note that the e-ph interaction constant G in 
Al is an order of magnitude larger than it is in noble 
metals (for comparison, Gai = 3.1 x 10 17 W/Km 3 vs. 
Gau = 2 x 10 16 W/Km 3 ), so that the energy exchange 
between the electron and phonon subsystems takes place 
on sub-picosecond times 18 . The rapid increase in the lat- 
tice temperature Tl and hence in the e-ph relaxation rate 
accounts for the fast rise of the pump-probe signal until 
the equilibrium between electron and phonon subsystems 
is reached at ~ 1 ps. This fast dynamics is superimposed 
with a slow decay due to the heat transfer from metal to 
glass substrate resulting in a plateau at 1-3 ps (see Figs. 
[2][3|). Importantly, the merging of electron thcrmaliza- 
tion and electron to lattice cooling stages implies that the 
short time dynamics is intrinsically nonequilibrium; it is 
therefore described below within the Boltzmann equation 
approach. 



The observed PM blueshift also indicates that the EOT 
is sensitive to the imaginary part of permittivity, e'^u)), 
for the following reason. Within the Drude model for 
Al, we have e' m (uS) = 1 — LofJ(uj 2 + 7 2 ). With such e' m , 
the SPP Bragg condition Eq. fl} leads to the following 
expression for, e.g., (0,1) mode frequency (at = 1): 
w sp P ~ t 1 - ( w 9 +7 2 )/ 2w p], where lu 9 = eg is the 
frequency associated with the reciprocal lattice period. 
Then an increase in 7 upon pulsed excitation leads to 
a weak redshift of the resonance frequency, in disagree- 
ment with the data. We thus conclude that the observed 
spectral blueshift cannot be explained by the present ver- 
sion of the SPP mechanism for the EOT phenomenon 2 . 
At the the same time, an increase in 7 causes an in- 
crease in the attenuation, e^, which for 7^0; can be 
approximated by ~ juj 2 /oj 3 . If e" n participates in de- 
termining the position and lineshape of the EOT bands, 
then the increase in 7 implies higher resonant frequen- 
cies (blueshift) in agreement with the data. The linear 
transmission is modeled below using effective medium ap- 
proach. 
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FIG. 4: (a) Transient spectra at different time delays are 
shown together with numerical fit (solid lines) based on solu- 
tion of the Boltzman equation in Sec. IHI Bl (b) Log derivative 
of measured zero-order transmission. 



III. THEORY 

The thansient optical spectra in metal films can be 
described in terms of time-dependent dielectric function 
e m (uj,t) — e m (uj) + 5e m (u),t), where e m (ui,t) is PM con- 
tribution due to a small (compared to ground state) frac- 
tion of photoexcited carriers. The differential transmis- 
sion AT(t) — T(t) — To can be then caculated from the 
linear transmission T (w) with s m (to, t) 14 ' 15 . While the 
time-resolved signal is primarily governed by the electron 
dynamics in metal film encoded in <5e m (w,i), the tran- 
sient spectra are determined by the effect of the latter on 
the dispersion of surface electromagnetic modes. Model- 
ing of the dynamics thus involves two distinct tasks: (i) 
description of the linear transmision through nanoholc 
array for real metal with absorptive dielectric function, 
and (ii) simulation of non-equilibrium electron-phonon 
dynamics in Al film. 

A. Effective medium model for linear transmission 

We begin with analytical model for the main transmis- 
sion band that incorporates the effect of the absorptive 
part of metal dielectric function. We are interested in 
the transmission of p-polarized light normally incident 
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on a metal film perforated with square array of cylin- 
drical nanoholes. We assume that the scatering in the 
direction normal to the plane of incidence is weak and 
adopt one-dimentional grating approximation for trans- 
verse magenetic and electric fields, H y (x, z) and E x (x, z). 
Magnetic fields in air (media 1), perforated metal (media 
2), and glass (media 3), can be written as 



H 



'Ml 



E-f ik 3n (z — h)+ingx 
^3n e , 



H 



2y 



t 2i e KlZ + r 2l e 



-Ki (h—z) 



(2) 
(3) 
(4) 



where kj n = ^Jk^tj — n 2 g 2 and kj = kjo (here fco = 
w/c), and rj n , and tj n are reflection and transmis- 
sion amplitudes of order n in medium j respectively 
(j = 1,3). In air/glass, the fields are expanded in the 
Bloch functions basis. In the medium 2 (perforated 
metal), the expansion goes over the normalized eigen- 
states of Maxwell equations (ME) in periodic medium, 
Hi(x)e~ KiZ , where eigenvalue Ki is the wave-vector along 
the normal to the surface. The eigenfunctions Hi(x) are 
localized in the hole region and are periodic over the 
array. The transverse electric fields in all three media 
are expressed through transverse magnetic fields via ME 
E x = —(i/koe)dH y /dz. In the relevant spectral range, 
it is sufficient to restrict to the lowest modes in air/glass 
(ri = 0,±1), and in perforated metal i = 1. 

The transmission coefficient is obtained by matching 
transverse electric and magnetic fields at air/metal and 
metal/glass interfaces. It is convenient to perform this 
matching in the Bloch basis for electric field, and in the 
basis of eigenstates Hi for magnetic field. This choice 
takes into account the fact that, for low volume frac- 
tion of holes, v = nd 2 /Aa 2 <C 1, the electric field is 
mostly reflected back into the air at the air/metal in- 
terface and transmitted to the glass at the metal/glass 
interface, while the magnetic field, conversely, mostly 
penetrates into the metal at the air/metal interface, and 
gets reflected back into the metal from the metal/glass 
boundary. We define magnetic and electric field moments 
averaged over the Brilloin zone (BZ) as 



(j nax Hi) BZ , 



H/BZ 



1 / ..v,^^, Hi (,x) 



e 2 (x) 



(5) 



BZ 



where £2(2) = f m + (th — £m)0H(x) is periodic dielec- 
tric function, e m (e/,.) being dielectric function in metal 
(hole) region, and Oh is the periodic step-function van- 
ishing outside the holes. Because we consider only the 
fundamental mode, we drop index i hereafter. By per- 
forming matching procedure at the interfaces z = and 



z = h we obtain the following system of equations 
t 2 +r 2 u = h (l + r w ) + 2hirn 

kin 



e„ (*2 - r 2 u) = 
t 2 u + r 2 = 

ks-n 

e n (t 2 u-r 2 ) = —^U n , 



, (5 n - r ln), 
^0^30 + 2/li^3i, 



with u 



£3^0 

z , yielding the zero-order transmission 
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Tn = \Un\ 2 = 



^12^236 



1 + ri 2 r 2 3e 



-2nh 



(6) 



(7) 



where t\ 2 = 2h 2 /(hi + n 2 ), £23 = 2n 3 /(h 2 + h 3 ) are ef- 
fective transmission amplitudes through the air/metal, 
and metal/glass interfaces respectively, and t\ 2 = (n 2 — 
n\)/{ni + n 2 ), r 23 — (n 3 — n 2 )/(h 2 + n 3 ) are the cor- 
responding effective reflection amplitudes. The effective 
refraction indices n are given by 



71 2 



nj\ 1 



2h\e\k 
h-oeok 



(8) 



where rij are refraction indices in air /glass (j — 1,3). 
The transmission coefficient is thus expressed via the 
overlap integrals, h n and e ra , between eigenstates in the 
periodic medium and in the air/glass. Such overlap 
causes hybridization of the states in perforated metal via 
the EM waves at the interfaces (both propagating and 
evanescent). It is this hybridization that determines the 
position and the lineshape of the transmission peak. 

In order to proceed further, we need to evaluate the 
overlap integrals, h n and e n , and the eigenvalue k. 
The former can be determined from the following self- 
consistency argument. While in the bulk of perforated 
metal the states Hi are strongly localized within holes, 
near the interfaces (where the overlap integrals are cal- 
culated) these states become delocalized due to the hy- 
bridization via EM in the dielectric. In other words, in 
the close vicinity of the interface, the profile of trans- 
verse magnetic field (that is continuous together with 
its derivative accross the interfaces) is similar at the 
both sides of the interface and can be approximated as 
H (x) = ho + 2h± cos gx. The relation between ho and hi 
on the metal side can therefore be determined from ME 
for H(x), 



d ( 1 8H\ 



)H = 0, 



(9) 



with periodically modulated dielectric function, e 2 (x) = 
60 + 2ei cos gx, where eo = e m (l — v) + ehf is the average 
and ei = v{t h - e m )£i, with £ n = sm(ngd/2)/(ngd/2), is 
the modulation amplitude 19 . In the first order in modu- 
lation, we obtain 



hi 
ho 



(10) 
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where ho is determined from the normalization condition. 
The moments of electric field, e n , are straightforwardly 
calculated from Eq. (0 as 

e = in(h e + 2h 1 e 1 )/k 

ei = in[h ei + hi(e + e 2 )]/k 0l (11) 

where l n = (e in <> x /e 2 (x)) B z = {e^ - f m >(„. 

Let us now turn to evaluation of the wavevector k that 
represents the lowest eigenvalue of ME Eq. ((9|. The 
eigenvalues are determined by the bulk states of the op- 
tically thick perforated metal film, and therefore are only 
weakly affected by the interfaces. In the bulk, the eigen- 
states are strongly localized within the holes and the 
weakly-modulated medium approximation does not ap- 
ply. In this case, the eigenvalues can be evaluated by us- 
ing the analogy between ME and Schrodinger equation 
(SE) with complex potential. Namely, for each region, 
Eq. © can be presented in the form d 2 H/dx 2 + (E — 
V)H = 0, where 

E = n 2 + k 2 [e{uu) - e h ] , 

V{x) = k 2 J2 [eH - e h ] 0(d/2 - \x + na\) , (12) 

n 

are the energy and potential, respectively. The latter is 
Kronig-Penney type attractive potential of periodically 
spaced quantum wells. Note that E = —k^Ch corresponds 
to k 2 — — £jn^o that is the value of k for a solid metal film. 
The fundametal mode is associated to the ground state of 
Eq. (|12[) with the energy spectrum being a band with the 
width determined by the tunnel transparency between 
the wells. In the case of high metal/dielectric contrast 
of e, so that ko<iy/\eh — e m \ 3* 1, this transparency is 
exponentially small and the ground state energy can be 
approximated by that of the lowest level in a single well. 
The latter can be found by replacing rectangular well 
by the delta- function potential, V(x) — Vq6(x), with the 
strength Vq — v(e m — eh)k 2 a, where we averaged over 
BZ. We then obtain E = —V 2 /4, yielding n 2 — —e m k 2 — 
u 2 a 2 k^{e m - e,0 2 /4. 

We thus arrive at the following physical picture for the 
main transmission band. Diffraction on an individual 
hole can be viewed as scattering on the bound state in 
the associated SE. The interference of diffracted waves, 
caused by the hybridization of localized hole states via 
EM waves in dielectrics, leads to the single-period mod- 
ulation of magnetic field profile at the interfaces. The 
latter is determined by the effective medium character- 
ized by weakly- modulated dielectric function depending 
on the hole volume fraction v. Note that eigenvalue k as 
well as magnetic/electric field moments h n ,e n are com- 
plex due to absorptive part of metal dielectric function 

In Fig. [T] we show the fit of linear trasmission obtained 
using the Drude form of dielectric function of Al, e m , with 
bulk plasmon frequency w p — 15.0 eV and damping rate 
7 = 0.2 eV— . It can be see that the main transmission 



band is described accurately, including the position of 
Wood's anomaly at u> s» 2.6 eV. Higher order maxima, 
including possible SPP peak at 4.0 eV are not captured 
within this model. 



B. Ultrafast Dynamics 

We now turn to the modeling of ultrafast dynamics. 
The transient spectra in a perforated film can be ob- 
tained from the linear transmission Tq(uj) by the replace- 
ment E m (u>) — > e m (ui) + S£ m (u>,t), where 5e m (u>, t) is 
the non-equilibrium contribution to the metal dielectric 
function. As noted above, 5e m (u),t) is mainly due to 
transient change in the e-ph relaxation rate "/(t) in the 
Drude dielectric function. The distinctive feature of elec- 
tron dynamics in Al films is a very fast (hundreds of fs) 
energy exchange between electron and phonon subsys- 
tems caused by strong e-ph interaction in Al. This en- 
ergy exchange takes place on the same timescale as the 
electron gas thermalization due to e-e scattering, render- 
ing the transient dynamics inherently non-equilibrium. 
Note however that the high-energy fraction of of phoex- 
cited electrons with the excess energy (E — Ep) ~ tko, 
scatter down on the timescale of tens of fs due to the 
quadratic energy dependence of e-e scattering rate, 7 ee <x 
(E — Ep) 2 , whereas the thermalization of electron pop- 
ulation is reached on the longer timescale of hundreds of 
f*^. Therefore, the higher energy end of electron dis- 
tribution relaxes primarily via e-e scattering processes 
while the energy transfer to the phonons primarily goes 
through the lower energy electrons. Since the overall 
heating of the lattice is of the order of Debye temper- 
ature (G>£> « 380K in Al), the phonons that gain energy 
from athermal electron gas can be considered as quasi- 
cquilibrium. The energy exchange between phonon and 
electron subsystems can be then described in terms of 
time-dependent latice temperature Ti,(t) as 

where Cl is the lattice heat capacity (Cx = 900 J/Kkg 
for Al), g(E) = E 1 / 2 (2m e ) 3 / 2 /2n 2 h 3 is the electron den- 
sity of states (DOS), and f(E,i) is the electron distribu- 
tion function that satisfies the Boltzmann equation 

Here the first two terms in r.h.s. are due to e-e and e- 
ph scattering processes respectively, and H(E,t) stands 
for electron-gas excitation by the pump pulse. Naturally, 
only the e-ph scattering term contributes to the energy 
exchange rate Eq. (|13p . Since the typical phonon energy 
is much smaller than the electron energy, the e-ph term 
can be more conveniently presented in the differential 
operator form^S. In the deformation potential approxi- 
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mation, the e-ph scattering term reads 



at 



i—ph 



Cep d 



f(l-f) + k B T L (t) 



df_ 

dE 



> (15) 



where C ep — tt 2 H 3 G / '^plhsffij and G is e-ph interaction 
constant. By substituting Eq. (fTo]) into (flU)) we obtain 
the rate equation for the lattice temperature, 



■ at 



where 



T;{t) 



G{T*-T L ) + S GL (t), (16) 



dEf(E,t)[l- f(E,t)] (17) 



plays the role of effective temperature of athermal elec- 
trons. Note that for thermal electron distribution /o, 
T* coinsides with the "regular" electron temperature T e . 
The last term in Eq. (fl7o|) is added to describe heat ex- 
change between the metal film and the glass substrate. 
With boundary conditions that the temperature and heat 
flux arc continuous accross the metal/glass interface, it 
has the formal: 5gl(*) — KcdTa/dz, where Tq and kg 
are the glass temperature and thermal conductivity, re- 
spectively {kq — 1-37 J/mK for SiC^). The above system 
should be supplemented by equation for heat transport 
in the glass, 



C G 



dT G 

at 



KG' 



3 2 T G 

dz 2 



(18) 



where Cq is the glass heat capacity (Cq = 1.63 x 10 7 
J/m 2 K for SiC>2). Note, however, that metal/glass dy- 
namics is important only for long times (see below). 

In the case short-time dynamics, the evolution of lat- 
tice temperature T^it) is governed by the interaction 
with strongly athermal electrons gas, and thus requires 
solution of Boltzmann equation for e-ph system (|14M16|) . 
We use the integral form for the e-e scattering term for 
parabolic band^S, 



aj_ 
at 



Es /2 EV 2 
Jx 



J dE l J dE 2 F(E u E 2 ,E 3 ) 



x 



1 + x 



arctan y/x 



(19) 



where K ee = rn e e 4 / (4h 3 e 2 ), e is static dielectric con- 

1/3 

stant of metal, and E s = (e 2 /27re)(3n e /7r) is the 
Thomas-Fermi energy. The phase space for scattering is 
determined by F(E 1: E 2 ,E 3 ) = f{E)f(E 1 )f(E 2 )f(E 3 ) - 
f(E)f(E 1 )f(E 2 )f(E 3 ), where f(E) = 1 - f(E) and 
E 3 = E + Ei — E 2 , while the angular integrals of 
Coulomb interaction yield the expression in the brackets 
in the integrand which is evaluated between the limints: 



u min,max 



max, mm 



(\/£iT^) ,(VEtVE^) IE, 



The last term in Eq. (114p . describing the energy deposited 
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FIG. 5: Non-equilibrium contribution to the electron distri- 
bution function in Al film is shown for several time delays. 
Calculations are based on solution of the Boltzman equation 
for the electron distribution function in Sec. IIII Bl 



by the pump pulse, is given by H (E, t) = Q(E)S(t) with 
Gaussian pulse profile 



S(t) 



(20) 



where Eq is the absorbed pump pulse energy density and 
(7l = 100 fs is the pulsewidth. The factor Q(E), describ- 
ing the phase space of phonon- assisted electronic transi- 
tions after absorption of a photon with energy u>, is avail- 
able in the literature 2 ^ and is included in the numerical 
calculations. 

In Fig. [5j we plot the time-dependent contribution to 
the electron distribution function, Af(E,t) = f(E,t) — 
f (E) , versus energy (measured from the Fermi level) at 
several time delays. For short times, electron distribu- 
tion is strongly non-equilibrium, as indicated by the long 
high-energy tail of A/. At t ~ 1 ps, electrons are almost 
thermalized and their distribution is close to the Fermi 
function of a hot electron gas. Note that, at 1 ps, the 
magnitude of Af has decreased significantly pointing to 
the energy transfer from electron gas to the lattice. 



IV. DISCUSSION 

The transient change in the lattice temperature, Ti,(i), 
provides the main contribution to the time-modulation of 
metal dielectric function, 



lu(uj + ijo) 1 



h(T L (t))- l0 ], (21) 



where j(u>,Tl) is the optical relaxation rate of an 
electron-hole pair excited by the probe light 22 (70 stands 
for the room-temperature relaxation rate). In the case 
when the lattice temperature remains much smaller than 
u>, the relaxation rate is nearly w-independent^ 7 -, and is 
given by 



7 (T) = B 



(22) 



7 



where J(y) = | dxx 4 coth^, and B is a material- 
dependent constant. In the relevant temperature range 
Tl ~ 6d, the Tt -dependence of 7 is nearly linear, 
i(Tl) ~ 7o + "f — Jo). Note that 70 contains contribu- 
tions also from the electron-impurity and e-e scattering 
processes; however those are not significantly affected by 
the pump pulse due to a relatively small change in the 
scattering phase space by the fraction of photoexcited 
electrons. 

The results of numerical simulations of short-time dy- 
namics are shown in Fig. [51 The calculated differential 
transmission dT/T is in excellent agreement with the 
data, reproducing the signal rise at subpicosecond times. 
In Fig. [6] we show the evolution of the lattice tempera- 
ture and the effective electron temperature. The rapid 
drop of T* after reaching its peak value at ~ 200 fs is ac- 
companied by increase in Tl during the period t < 1 ps, 
indicating energy intake by the phonon subsystem while 
the electron gas is still in the thermalization stage. After 
electron and phonon subsystems equilibrate, the signal 
changes weakly (plateau) until the slower process of en- 
ergy transfer to the glass substrate takes over, leading to 
the decay of pump-probe signal from the metal [see Fig. 
[3j. For long time delays, electron and phonon subsystems 
are in thermal equilibrium with each other, and the long- 
time dynamics is well descrbed by the two-temperature 
model that incorporates the heat exchange between the 
metal and the glass^i. To compare to Boltzmann equa- 
tion approach, here we also included in calculations the 
rate equation for electron temperature T e (t) that can be 
straightforwardly obtained from Eqs. (|I4M1 6[) by the re- 
placement of f(E,t) with thermal distribution: 

C e 8T e /dt = G(T L - T e ) + S{t), (23) 

where C e = (3T e (with f3 = 134 J m" 3 KT 2 for Al) is 
the electronic heat capacity. After ~ 40 ps, the signal 
reaches a plateau corresponding to thermal equilibrium 
between film and glass [see Fig. [3J. Here the evolution 
of dT/T also follows mainly the phonon dynamics [see 
Fig. [f^b)]. Note, however, that the subpicosecond sig- 
nal rise is noticably faster than the one calculated using 
Boltzmann equation approach [cf. Fig. [5]. The reason 
is that thermal distribution overestimates the fraction of 
lower energy electrons and, hence, the rate of their energy 
transfer to the lattice. The measured signal also exhibits 
weak oscillations with a period ~ 5 ps; those could be 
attributed to the vibrational mode accross the film that 
is excited by the rapid heating of the lattice, but are not 
captured by the present model. 

The calculated transient spectra are shown in Fig.[5]for 
several time delays in the long-time regime when the elec- 
tron and phonon subsystems are in thermal equilibrium 
with each other. The blucshift is well reproduced as well 
as is the nonlinear cj-dependence of dT/T in the main 
transmission band frequency range. The origin of the 
blueshift can be understood as follows. The eigenmodes 
in a nanohole array are the result of interference of multi- 
ple waves scattered from periodic metal/dielectric regions 




0.0 0.5 1.0 1.5 
Time delay (ps) 




3.2 - 

_J , , , , , , , , , 

10 20 30 40 

Time delay (ps) 

FIG. 6: (a) Short-time evolution of effective electron (solid 
line) and lattice (dashed line) temperatures in Al film fol- 
lowing the excitation, (b) Long-time evolution of the lattice 
temperature due to heat transfer to the glass matrix. 



as described by the above effective medium model. In the 
presence of attenuation, the phase of scattered wave has a 
contribution cx Re-y/e^ = e" n / \/\s' m \ w jUp/ui 2 . There- 
fore, a transient change in attenuation causes correspond- 
ing shifts in scattering phases. This, in turn, changes the 
interference of waves making up an eigenmode. To com- 
pensate for the phase shifts, a transient increase in 7 
must be accompanied by a corresponding increase in u>; 
and consequently the transmission resonances blueshift. 



V. CONCLUSION 

In summary, we investigated the ultrafast dynamics 
of Al-based plasmonic lattice that reveals a non-trivial 
role played by the plasma attenuation in determining 
the EOT resonant bands. While its effect on the SPP 
dispersion is relatively weak, attenuation is more impor- 
tant for modes comprised of a multitude of scattered 
waves. By directly altering their interference, ultrafast 
time-resolved optical techniques unravel the unique prop- 
erties of photoexcited SEW in plasmonic metals. 

This work was supported in part by NSF DMR Grant 
No. 05-03172, by NSF Grant No. DMR-0606509, and by 
NIH Grant No. 2 S06 GM008047-33. 
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